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Abstract 

We investigate the CP violating asymmetry, the forward backward asymmetry and the CP 
violating asymmetry in the forward-backward asymmetry for the radiative dileptonic B-meson 
decays B& — > 'y£ + £~for the £ = e, /i, r channels. It is observed that these asymmetries are 
quite sizable and — ► 7 decays seem promising for investigating CP violation. 



1 Introduction 

The rare B meson decays are one of the most important probes of the effective Hamiltonian 
governing the flavor-changing neutral current transitions b — > s(d)£ + £~. Among them, the rare 
Bs,d — ► 7 £ + £~ decays receive a special interest due to their relative cleanliness and sensitivity to 
new physics. They have been investigated in the framework of the SM for light and heavy lepton 
modes in refs.[JT]]-[|5[]. These decays have also been studied in models beyond the SM, such as 
MSSM [§] and different versions of the two Higgs doublet models (2HDM) and it was 

reported that the new physics effects can give sizable contributions to the relevant observables. 

For b — > s£ + £~ transition, the matrix element contains the terms that receive contributions 
from ti, cc and uu loops, which are proportional to the combination of £j = V t bV t *, £ c = V^V^,* 
and £ u = V u bV* s , respectively. Smallness of £ u in comparison with £ c and £ t , together with 
the unitarity of the CKM matrix elements, bring about the consequence that matrix element for 
the b — > s£ + £~ decay involves only one independent CKM factor £ t , so that the CP violation 



in this channel is suppressed in the SM [11, 12]. However, for b — > d£ + £ decay, all the CKM 



factors rjt = V t bV t * d , tj c = V c t>V* d and rj u = V u bV* d are at the same order in the SM so that they 
can induce a CP violating asymmetry between the decay rates of the reactions b — > d£ + £~ and 



b — > d£ + £~ [13]. So, b — > d£ + £~ decay seems to be suitable for establishing CP violation in B 
mesons. On the other hand, it should be noted that the detection of the b — ► d£ + £~ decay will 
probably be more difficult in the presence of a much stronger decay b — ► s£ + £~ and this would 
make the corresponding exclusive decay channels more preferable in search of CP violation. In 
this context, the exclusive B^ — ► p£ + £~ and B^ — > 7r£ + ^~ decays have been extensively studied 
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in the SM [14] and beyond [ 15]-[jT^]. So, we think that it would be interesting and complementary 
to consider the remaining exclusive mode Bd — > 7 £ + £~ . 

In this paper, we would like to study the CP violation in the exclusive B d — > 7 ^ + ^~decay in 
the context of the SM. B& — > 7^ + ^~decay is induced by the pure-leptonic decay Bd — > ^ + ^~, 
which is well known to have helicity suppression for light lepton modes, having branching ratios 
(BR) of the order of 10 -15 for £ = e and 10~ 10 for I = // channels [Eh. However, when a 
photon line is attached to any of the charged lines in Bd — > process, it changes into the 

corresponding radiative ones, Bd — ► "yl + £~, so helicity suppression is overcome and larger 
branching ratios are expected. In [||]), it was found that in the SM, BR{Bd — ► £ + ^~7) = 
(1.5(1.5) , 1.2(1.8) , - (6.2)) x 1CT 10 for £ = e,/x,r, respectively. Although these BR's are 
quite low, in models beyond the SM they can be enhanced by two (one) orders, as shown e.g. in 
[20]([[2~I|]) for B s / d ) ~> 1 ^ £~ decay, so investigation of this process may also be interesting from 



the point of view of the new physics effects. 

In Bd — > 7 £ + ^~decays, depending on whether the photon is released from the initial quark or 
final lepton lines, there exist two different types of contributions, namely the so-called "structure 
dependent" (SD) and the "internal Bremsstrahlung" (IB) respectively, while contributions coming 
from the release of the free photon from any charged internal line will be suppressed by a factor 
of ml/Myy- The SD contribution is governed by the vector and axial vector form factors and it 
is free from the helicity suppression. Therefore, it could enhance the decay rates of the radiative 
processes Bd — ► l + £~^ in comparison to the decay rates of the pure leptonic ones Bd — ► £ + £~ . 
As for the IB part of the contribution, it is proportional to the ratio me/ tub and therefore it is 
still helicity suppressed for the light charged lepton modes while it is expected to enhance the 
amplitude considerably for £ = r mode. However, we note that IB part of the amplitude does 
not contribute to CP violating asymmetry Ac p and the forward-backward asymmetry App (see 
section 2). 

We organized the paper as follows: In section ||, first the effective Hamiltonian is presented and 
the form factors are defined. Then, the basic formulas of the differential branching ratio dBR/dx, 
Acp, Apb and CP violating asymmetry in forward-backward asymmetry Acp(App) for Bd — > 
7 £ + £~ decay are introduced. Section || is devoted to the numerical analysis and discussion. 

2 The theoretical framework of Bd —> 7 ^"decays 

The leading order QCD corrected effective Hamiltonian which is induced by the corresponding 



quark level process b — > d£ + £ , is given by []22[] - [p5[] : 

7 1 Gu)[O?(a0-Oi0u)] 



AC ( 10 



+C 2 ( M )[0£( M )-0 2 (//)]} (1) 



where 



u ~ T/ V* ' y ' 

V tbV t d 

using the unitarity of the CKM matrix i.e. V t bV t * d + V ub V* d = -V cb V* d . The explicit forms of 
the operators Oi can be found in refs. [||, ||]. In Eq.@, d(fi) are the Wilson coefficients 



calculated at a renormalization point \i and their evolution from the higher scale /x = mw down 
to the low-energy scale \i = m b is described by the renormalization group equation. For Cj (n) 



2 



this calculation is performed in refs.[|(], 27] in next to leading order. The value of Cio(m&) to the 
leading logarithmic approximation can be found e.g. in [E2|, 25]. We here present the expression 



for Cg(/x) which contains the terms responsible for the CP violation in Bj — > 7£ + £~decay. It 
has a perturbative part and a part coming from long distance (LD) effects due to conversion of the 
real cc into lepton pair 



C e 9 ff ^) = Cr\ri+Y reson {s) 



(3) 



where 



sipert 
°9 



0) = C 9 + h{u,s)[ZC 1 { l i) + C2{^ + ?,C^) + C^)+ZC^) + C & {i J ) 
+ A u (3d + C 2 )\ - h(l, s) (4C 3 (/x) + 4C 4 (/x) + 3C 5 (m) + C 6 0)) 

- \h(0, s) [C3Q1) + 3C 4 (aO + A u (6Ci(At) + 2C 2 (/i))] 



(4) 



and 



Y 

1 r( 



x [(3d( M ) + C 2 (/i) + 3C 3 (m) + C 4 (/i) + 3C 5 (/i) + C 6 ( M )) 



(5) 



In Eq.(f|), s = q 2 /m 2 B where q is the momentum transfer, u = ^ and the functions h(u, s) arise 
from one loop contributions of the four-quark operators 0\ — Oq and are given by 



h(u, s) 



8, m b 8, 8 4 

-— In lnu H 1 — y 

9 /x 9 27 9 



h(0,s) 
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8 m b 
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In s H — Z7r . 
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for y 



< 1 



fory = ^>l, 



(6) 



(7) 



The phenomenological parameter k in Eq. (g) is taken as 2.3 (see e.g. [JT3[]). 

Neglecting the mass of the d quark, the effective short distance Hamiltonian for the b 
d£ + £~ decay in Eq.([j]) leads to the QCD corrected matrix element: 



M 



Gfol 



V tb V t *A C 9 e// (m 6 ) d 7M (l - 75)6*7"* + C 10 (m 6 ) d 7/1 (l - 75 )6 ^75^ 



2^2 



7T 



- 2C 7 e// (m 6 ) ^dia^l + 75)6*7^ 



(8) 



The next step is to calculate the matrix element of the B^ 
the sum of the SD and IB parts 



0+0- 



decay. It can be written as 



M(B d ^£+r-f) 



M S D + Mib ■ 



(9) 
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It is evident from Eq.^ that the following matrix elements are needed for the calculations of 
M S d part : 

< 7 |<i7 M (lT75)&|5) = A r {e^e*»q x k°g(q 2 )±i\e;(k q )-(e*q)k fl ]f(q 2 )}, 



(7\dia^q»(lTj 5 )b\B) = e^e^k* 9l (q 2 ) =p i\e*Jkq) - (e*q)kj fi(q 2 )W) 



m 



B 



where e* and are the four vector polarization and four momentum of the photon, respectively, 
and pb is the momentum of the B meson. The form factors g, f, g±, and f± were calculated in the 
framework of the light-cone QCD sum rules in [||, |], [28[ ]. 



M 



The mauix element describing the SD part can be obtained from Eqs. (|8|)-(|10|), 



SD 



. nr VtbVtd — 

4V2 it m z 



B 



£-f(l - l5 )l [Axe^e^q^ + iA 2 (e;(kq) - (e*q)k, 



+ l-fil + Ts)* \Bie^ a pe* v q a kP + iB 2 (e*Jkq) - (e*q)k 



(11) 



where 



A x = 


-4m b C e 7 ff gi + (C e 9 ff 
Q 


- C 10 )g , 


A 2 = 


=^m b C e 7 ff h + {C e 9 ff 
qz 


- c w )f , 


Bi = 


^-m b C e 7 ff gi + (C e 9 ff 
qz 


+ Cva)g , 


B 2 = 


=4m b C e 7 ff h + (C e 9 ff 


+ C w )f , 



For the IB part, using 

(0\d W b\B) : 

and conservation of the vector current, we get 

oGf 



SbVb^ , (0\da^(l+ l5 )b\B) = 



M 



IB 



Ay/2 7T 



VMef B i{F> 



J B 



2pik 2p 2 k 



75* 



with 



4m £ Ci 



(12) 



(13) 



(14) 



(15) 



Substituting Eqs.(ll ) and ( Jl4| ) into Eq. @, squaring it and then averaging over the initial and 
summing over the final spins of the leptons and polarization of the photon, we find the photon 
energy distribution, after integration over the phase space, which is given by 



dT 

dx 



aG F 2 
VtbVtd 



4^2 vr 



a 



(2tt)' 



irrriB D(x) 



(16) 
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where 



D(x) 



2 3 

m B x v 



I(|A 1 | 2 + |A 2 | 2 + |Bi| 2 + | J B 2 | 2 )(l-r-x) 



+r Re(A 1 B* + A 2 B*) 



+ f B m e x 2 Re([A 1 + B X \F*) In 



l + v 
1-v 



-I 



2v- 



1 - x 



+ 2 + 



4r 



x In 



x x 



l + v 
1-v 



(17) 



where x = 2E~ j /m B is the dimensionless photon energy and v = y 1 — with r = m 2 /m? B . 

We now consider the CP violating asymmetry, Acp, between the B^ — > ^f£ + frm& B^ — > 
7 £ + £~ decays, which is defined as follows: 



A CP {x) 



-)-r(B d - 7 *+r 



T(B d -> 7 «") + r(5 d -> 7 W 
Using this definition we calculate the >1cp as: 



f(D(x) - H(x)) dx 



(18) 



(19) 



where 



H(x) 



-2x 2 



Im A, 



2 im £ 2 I Cj ff (f fi+g gi)m h v x 



x — 2r — 1 
1 — x 



+6C w gfBm e In 



2, 1 + V 



1-U 



(/ 2 + g 2 )m B vx(x — 2r — 1) Im £^£2 



In calculating this expression, we use the following parametrization: 

c 9 e// = 6 + A M 6 . 



(20) 



(21) 



We note that in these integrals the Dalitz boundary for the dimensionless photon energy x is 
taken as 



5 < x < 1 



Ami 



m 



2 ' 



(22) 



since |-M/b| term has infrared singularity due to the emission of soft photon. In order to obtain 
a finite result, we follow the approach described in ref . |4j] and impose a cut on the photon energy, 
i.e., we require E y > 25 MeV, which corresponds to detecting only hard photons experimentally. 
This cut requires that E y > 5m B /2 with 5 = 0.01. 

Next, we consider the forward-backward asymmetry, Ap B , in B^ — > ~/l + £~. Using the 
definition of differential Ap B 



[ l r1z^- - f° dz^- 

JO az dz J -I az dz 

r x d z £- + f° dz^ ' 

JO aZ dz + J-l aZ dz 



Afb(x) 



(23) 



where z = cos 9, 9 is the angle between the momentum of the B-meson and that of I in the cm. 
frame of the dileptons £ + l~, we find 



A 



FB 



J dxE(x) 
j dxD{x) 



(24) 
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with 



E(x) = -Avx 2 (m% x sjjx - l)(x - 1 + Ar)Re(A 1 A* 2 - B\B 2 ) 

( x — 1 \ 4r \ 
+ 4f B m e v — — In -Re((A 2 - B 2 )F*) . 



1 + 4r 



(25) 



We have also a CP violating asymmetry in App, Ac p (Afb), which is an important measur- 
able quantity in extracting precise information about free parameters of the models used. Since 
in the limit of CP conservation, one expects Afb = —Afb I^L where Afb and Afb are the 
forward-backward asymmetries in the particle and antiparticle channels, respectively, it is defined 

as 

Acp(Afb) = A FB + A FB , (26) 
Here, Afb can be obtained by the replacement, 

C 9 e// (A U ) -> C e J f (\ u -» A;). (27) 



3 Numerical analysis and discussion 

In Figs. (JTJ-[5|), we present the dependence of the Aqp, Afb and Acp(Afb) on the dimensionles 
photon energy x for the — ► 7 £ + £~(£ = e, p, r) decays for two different sets of parameters 
(p, 77) = (—0.07; 0.34) and (0.3; 0.34) in the following Wolfenstein parametrization: 

_ p(l - p) - y 2 - iy , 2) 
A «" (l- p) 2 + r/ 2 +°( A )- (28) 

We have also evaluated the average values of CP asymmetry < Aqp >, forward-backward asym- 
metry < Afb > and CP asymmetry in the forward-backward asymmetry < Ac p (Afb) > in 
Bd ^ 7 ^ + ^~decay for the above sets of parameters (p, 77), and our results are displayed in Table 
1 and 2 with and without including the long distance effects, respectively. 

For the form factors g, f, g\ and f\, we have used the values calculated in the framework of 
light-cone QCD sum rules in refs. |], 28], which can be represented in the following dipole 
forms, 

^ . 9(0) ,,,, no) 



9(<f) = , : N2 , /(<r) — ^, 



(i-4Y 



/ 2n 51(0) , 2^ /l(°) , ocn 

) = 72 » A (« ) = 7T » ( 29 ) 

,2 

/1 



1-4- 1-4 

ml / V m 



where 



5 (0) = 1 GeV , f(0) = 0.8 GeF , g\ (0) = 3.74 GeV 2 , fi (0) =0.68 GeV' 
m g = 5.6 GeV , m f = 6.5 GeV , m gi = 6 AGeV , m fl = 5.5 GeV . 
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0; if) 


< A C p > 
£ = e I = ii I = t 


< A C p{A FB ) > 
£ = e £ = ix £ = t 


< A FB > 
£ = e I = [jl £ = t 


(0.3; 0.34) 


0.118 0.114 0.103 


-0.004 -0.004 0.038 


-0.472 -0.460 -0.188 


(-0.07; 0.34) 


0.061 0.059 0.050 


0.004 0.004 0.021 


-0.503 -0.490 -0.192 


Table 1: The average values of A C p, A FB and A C p{A FB ) in B d — > 
lepton modes including the long distance effects. 


7£ + £~for the three distinct 


(p; v) 


< A C p > 
£ = e £ = f± £ = t 


< A C p(A FB ) > 
I = e £ = f± I = t 


< A FB > 
£ = e £ = ji £ = t 




(0.3; 0.34) 


0.106 0.103 0.081 


0.011 0.010 0.021 


-0.520 -0.507 -0.200 




(-0.07; 0.34) 


0.053 0.051 0.034 


0.010 0.010 0.011 


-0.529 -0.514 -0.200 





Table 2: The same as Table (P, but without including the long distance effects. 



In addition to these form factors, the input parameters and the initial values of the Wilson coeffi- 
cients we used in our numerical analysis are as follows: 

m B = 5.28 GeV , m b = 4.8 GeV , m c = 1.4 GeV , f B = 0.2 GeV , 

m T = 1.78 GeV, m M = 0.105 GeV, \V tb V t * d \ = 0.01, 

Ci = -0.245, C 2 = 1.107, C 3 = 0.011, C 4 = -0.026, C 5 = 0.007, 

C 6 = -0.0314, Cj ff = -0.315, C 9 = 4.220, C w = -4.619. (30) 

In our numerical analysis, we take into account five possible resonances for the LD effects 
coming from the reaction b — > dtfi — » d£ + £~, where i = 1, 5 and divide the integration region 
into two parts for t = r: 5 < x < 1 - ((m^ 2 + 0.02)/m B ) 2 and 1 - ((m^ 2 - 0.02)/m B ) 2 < 
x < 1 — (2m^/ms) 2 , where m^ 2 = 3.686 GeV is the mass of the second resonance and into 
three parts for I = e and ii: 5 < x < 1 — ((m^ 2 + 0.02)/ms) 2 , 1 — ((m^ 2 — 0.02)/m#) 2 < 
x < 1 - ((m^ + 0.02)/m B ) 2 and 1 - ((m^ - 0.02)/m B ) 2 < x <l - (2ra f /ra B ) 2 , where 
m^, 1 = 3.097 GeV is the mass of the first resonance. 

For reference, we present our SM predictions without long distance effects 

BR(B d ^j£ + r) = (8.43,8.52,6.22) x 10" 10 , (31) 

for £ = e,fi,T. Here, we have used r(B d ) = 1.5 x 10~ 12 s, |V^V^| = 0.01, and (p;rj) = 
(0.30; 0.34). Our result for £ = r is in a good agreement with ref.[^], but for £ = e and /i, they 
are larger than those of ref.[Q] and 

In Fig.([l|) and Fig.(f|), we present the dependence of Aqp on the dimensionless photon energy 
x, for B d — » 7^ + £~decay for the Wolfenstein parameters (p; rf) = (—0.07; 0.34) and (p; rj) = 
(0.3; 0.34), respectively. The three distinct lepton modes £ = e, p, r are represented by the 
small dashed, dashed and solid curves, respectively. We observe that the Acp for £ = e, p cases 
almost coincide, reaching up to 15 % for the intermediate values of x. The Aqp for £ = r mode 
exceeds the values of the other modes and reaches 40 % in the high-x region. We also observe 
from Tables 1 and 2 that including the LD effects in calculating < Aqp > changes the results 
only 11 — 15% for £ = e, p modes, while £ = r mode, it is quite sizable, 27 — 47%, depending 
on the sets of the parameters used for (p; rj). 
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The x dependence of Afb for the B^ — > 7 £ + £~(£ = e, p, r) decays are plotted in Fig.([|) 
for (p; 77) = (0.3; 0.34). Since we observe that this dependence is almost unchanged for the 
other set of (p; 77) = (—0.07; 0.34), we do not display it here. We see that Afb is negative for all 
values of x, except in the resonance regions. < Afb > amounts to —50% for £ = e, fi modes, 
but it stands smaller for £ = r mode, (—20%) as expected. The LD effects on < Afb > are 
about 10%, but in reverse manner, decreasing its magnitude in comparison to the values without 
LD contributions. 

We present the dependence of the Acp(Afb) of — ► 7 decay on x in Fig.(||) and 
Fig.(|5|) again for two different sets of the Wolfenstein parameters. As for Aqp, Acp(Afb) for 
£ = e, and I = /j modes almost coincide. They can have both signs and stand smaller for all 
values of x compared to the one of £ = r mode. For the latter case, Ac p (Afb) is positive for 
all values of x except in the resonance region and it is at the order of magnitude 10 _1 . LD effects 
seem to be quite significant for < Ac p (Afb) >, decreasing its value by 60 — 80% for £ = e, p 
modes, while £ = r mode, including LD effects increases < Acp(Afb) > by 75 — 90%. 

As a conclusion we can say that there is a significant Acp and Ac p (Afb) for the B^ — > 
7£ + £~decay, although the branching ratios predicted for these channels are relatively small be- 
cause of CKM suppression. Experimentally, to be able to measure the BR of B^ — > 7 £ + £~decays, 
the required number of events, N, are N ~ BR x (number of B — mesons produced). Since 
approximately, 6 x 10 11 B^ mesons are expected to be produced per year, we may hope that 
Bd —>■ 7 £ + £~ decays could be measured in LHC-B experiment in future. 
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Figure 1: A C p for B d — > 7 £ + £~ decay for the Wolfenstein parameters (p, 77) = (—0.07; 0.34). The 
three distinct lepton modes £ = e, /i, r are represented by the small dashed, dashed and solid curves, 
respectively. 




Figure 2: The same as Fig. ([I]) but for the Wolfenstein parameters (p, rf) = (0.3; 0.34) 
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Figure 3: Afb for — > 7^ + £~decay for the Wolfenstein parameters (p, rf) = (0.3; 0.34). The 
three distinct lepton modes £ = e, p, r are represented by the small dashed, dashed and solid curves, 
respectively. 
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Figure 4: A CP (A FB ) for B d — > 7 £ + £~decay for the Wolfenstein parameters (p, rf) = (—0.07; 0.34). 
The three distinct lepton modes i — e, /i, r are represented by the small dashed, dashed and solid 
curves, respectively. 
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Figure 5: The same as Fig. ([I]) but for the Wolfenstein parameters (p, rj) = (0.3; 0.34) 
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